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Abstract 

We study the scalar perturbation sector of the general axisymmetric warped 
Salam-Sezgin model with codimension-2 branes. We focus on the perturbations which 
mix with the dilaton. We show that the scalar fluctuation analysis can be reduced 
to studying two scalar modes of constant wavefunction, plus modes of non-constant 
wavefunction which obey a single Schrodinger equation. Prom the obtained explicit 
solution of the scalar modes, we point out the importance of the non-constant modes 
in describing the four dimensional effective theory. This observation remains true 
for the unwarped case and was neglected in the relevant literature. Furthermore, we 
show that the warped solutions are free of instabilities. 



mail address: hyun.min.lee@desy.de 

mail address: antonios. papazoglou@epfl.ch 



1 Introduction 



Among six dimensional supergravities, the Salam-Sezgin model [1] (the supersymmetric 
analogue of [2]) has received particular attention for the past decades. It has the attractive 
feature that it gives rise to a massless chiral gravitino in four dimensions, thus reducing 
the supersymmetry to = 1. By enlarging the gauge group and adding a number of 
hypermultiplets [3] the model can be rendered anomaly free-*^ via the Green-Schwarz mech- 
anism. The later model is more relevant to phenomenology since it contains a large gauge 
group with matter fields charged under it. Recently, there has been renewed interest in 
the Salam-Sezgin model, due to the fact that the vacua that are obtained, exhibit the 
mechanism of spontaneous compactification with gauge field fluxes. 

Another particular characteristic of this model, that was noticed recently [5], is that 
all the non-singular {i.e., with no singularities worse than conical) maximally symmetric 
vacua are of the type (Minkowski)4 x X2, with X2 a two dimensional manifold. The 
generic vacuum solutions of this type have been found to get a warping in front of the four 
dimensional line element [5-7]. This warping leads always to the appearance of conical 
singularities which have to be supported by codimension-2 branes. The unique vacuum 
that preserves A^ = 1 supersymmetry is the one that has no warping (and with no branes 
present) in the four dimensional world-volume and that has its gauge field flux embedded 
in the gauged U{1)r direction. The appearance of the warping breaks the remaining 
N — 1 supersymmetry. More general warped solutions without axial symmetry of extra 
dimensions have been found in Ref . [8] . 

The importance of such kind of compactifications has been increased the last years with 
the consideration of models which try to ameliorate the cosmological constant problem, 
the so called selftuning models [9]. In these models, use was made of a special property of 
the codimension-2 branes, that they do not curve their world-volume, but instead induce 
a deficit angle in the bulk [10]. Thus, the vacuum energy of fields hving on these branes 
does not gravitate in four dimensions. This would give the hope to solve the puzzle of 
the smallness of the cosmological constant. However, from the available models, the ones 
with flux compactiflcations have a hidden flne-tuning related to the flux quantization or 
conservation condition [5, 7, 11, 12]. The cosmology on a thick codimension-2 brane in six 
dimensional flux models has been studied in Ref. [13] and has come to the same conclusion 
regarding the fine-tuning. Recently, different kinds of compactifications with de Sitter 
(singular) 3-branes [14] or a 4-brane only [15] have been put forward in the Salam-Sezgin 
model. 

In the present paper, we will discuss the linearized scalar fiuctuations of the Salam- 
Sezgin model for the general warped background of the form given in Refs. [5,7]. A similar 
analysis for a non-sup ersymmetric model has been recently done in [16]. Although we focus 
on the warped solution with axial symmetry of extra dimensions, our fluctuation analysis is 
also applied to the general warped background without axial symmetry in the local patch 
coordinate for each brane. We will study only the fiuctuations which are coupled with the 

^See [4] for recent advances on constructing anomaly free models in six-dimensional supergravities. 
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dilaton perturbations. We see that they are divided to fluctuations with constant profile 
along the extra dimensions and to fluctuations with non-trivial wavcfTinctions. The lowest 
massive non-constant mode mixes however with the massive constant mode. This mixing 
is always present irrespective of the presence of warping and has been neglected in the 
hterature when discussing the effective four-dimensional physics of the model [5, 12, 17]. 
In particular, in the four-dimensional supergravity description of the unwarped solution 
without branes [17], the new non-constant mode has to be included as a new massive 
chiral multiplet relevant for low energy physics. Moreover, we can see that the solutions 
of general warping are free of instabilities. 

in the warped case, we have found the interesting result that the mixing with the new 
mode plays a crucial role in determining the instability of the solution, from the wrong 
sign of the kinetic term of one mode for some region of the parameter space. 

The paper is organized as follows. In the beginning we will review the Salam-Sezgin 
model and its vacua in general (warped and unwarped) . In the following we will derive the 
perturbation equations and reduce them to a single Schrodinger. Then, we will present 
the explicit wavefunctions and masses for scalar modes. Next, we will derive the effective 
action, showing the mixing of the constant massive mode and lowest massive non-constant 
mode and comment on the stability of the warped solution. Finally, the conclusions will 
be drawn. The detailed derivation of the linearized equations and the quadratic effective 
action for the scalar modes is presented in the following two appendices. 



2 Salam-Sezgin model: general axisymmetric vacua 

We will first review the general vacuum solutions of Salam-Sezgin model [1] with axial 
symmetry. The bosonic sector of the system consists of the metric Qmn^ a Kalb-Ramond 
field Bmn, a dilaton and a gauge field Am- For the purpose of this paper, we will set the 
Kalb-Ramond field to its zero background value, so we will not include it in the action. 
Then, the bosonic bulk action of the system is given by 



— fiq e 2'. 



supplemented with the 3-brane action as 

'S'brane j d'^X^/^Vs = j d^X^J^L/^ , 

where a distributional brane energy density is 



-9^ 



U = - \ d''x^l-^Vs6^''\X -X{x)) 



(2) 



(3) 



Here, Vg is the tension, is the metric pulled back to the brane worldvolume and X^(x) is 
the embedding of the brane in the six-dimensional bulk. The gauge coupling g corresponds 
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to the gauged U{1)r of the model and in principle different from the the gauge coupling g 
of the gauge field Am- The Einstein and field equations derived from the above action are 



Rmn = e-^'f'gMN + ^e'^'^iFupFN ^ - ^9mnF^q) 

+ \dMCl>dN(l> + f'MN : (4) 

with T^N being the brane contribution and 

lei^F|Q-8/e-i^ (5) 
dMiV^e'^'^F^'') = . (6) 

Assuming the axial symmetry in the internal space, there will be in general, two 3-branes 
sitting in the antipodal points of the axis of symmetry. The general warped solution in 
this case can be analytically found in the following gauge [5, 7] 



ds^ = W^{r)rj^,dx^'dx'' + ^\r)[dr^ + X''a''{r)d9^] , (7) 
(/) = UnW , (8) 

with the various functions given by 

JO Ji 
W - fy, /o = l + ^, /l = l + ^, (11) 



n 



where flux is a constant and the two radii are given by 

The quantization condition of the gauge field fiux is given by the relation 

— n, n = mteger . [16) 

For the given metric solution, the brane contribution to the Einstein equation is ex- 
pressed by 
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fl,^ = ~^Vs (g,. d'i.d^^ - gMN)S^'\y- Vs) , (14) 
^ V 9 

where ^5 is the brane position. The definition of the delta function is such that J d^y5'^'^\y) = 
1, and thus if expressed in the above coordinates (with r the radial coordinate), it is 
S^'^Kv) = ^(^)/(2'7r). In this general solution, the metric has two conical singularities, one 
at r = and the other at r = oo^, with deficit angles Ss (supported by tensions Vg — 2Ss) 
given by 



2n 



1-A, 
1 



^2 



(15) 
(16) 



For tq — ri, i.e., for q — Ag, we have the unwarped model. In the case when A = 1 
and g — g, the unwarped case is possible only if n = 1, i.e., with no branes present. In all 
cases where 3-branes are present in the vacuum, supersymmetry is completely broken. 

Finally, let us go to a coordinate system which is Gaussian-normal with respect to the 
two branes. In this new radial coordinate the perturbation equations in the next section 
will be expressed in a more convenient way. Thus, if we define 



dp — ^ydr, a — , 



the metric is expressed as 



dsl = W'^^dx^dx"" + dp' + X^a^de 



(17) 



(18) 



The background equations of motion are expressed in the new coordinate p as (all 
primes from now on are derivatives with respect to p) 



W 



12 



+ 



W'a' 

Wa 
W'a' 

Wa 

W'a' 



"a ^ Wa 



where we have substituted 



W" 

+ 



w 
1 

w^ 



2 Q 

-9 + 



1 



2/ + 



3q^ 
8W^ 



(j)^4:\nW and Fp, 



AirXa 



Kp - Ps 



(19) 
(20) 

(21) 



\qa 
W^ 



(22) 



Note that this point is at a finite proper distance from r = 0. 
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In the unwarped case we can make this change of coordinate exphcitly and have that 

a(rt^^.n(!^). (23) 

3 Linearized scalar perturbations 

We would hke in this section to perturb the above general vacuum solution and in particular 
the spin-0 sector. For this purpose we consider the following ansatz for the perturbed 
metric: 

dsl = e-'^W\,dx''dx'' + e^(dp2 + ^2(^-i) x'^a'^dO^) . (24) 

The perturbation in front of d6^ is the right one to avoid mixing with the graviton - or 
in other words it comes from the (pLu) Einstein equations with n ^ v (see e.g., [18]). The 
gauge field perturbation is considered as following 

F^,e = V^ae , Fpe = ^ + a!^ , (25) 
with all the other components vanishing. In addition, the scalar field perturbation is 

4) = AhiW + f. (26) 

In all the above perturbations ip, ^, /, ae are functions of the 4d coordinates {x) and 
the radial one (p). We will not consider the ^-dependence which will provide the angular 
excitations of the resulting modes, i.e., we will restrict ourselves to the s-mode of the 
excitations. 

The above scalar perturbations are a subset of the most general perturbations. The 
most general perturbations would include also Ap, A (with — d/^A), = Qpo, b = Bpg 
and B (the scalar dual of -B^jy) in addition. However, as it was shown in [19], the complete 
set of perturbations is divided into two subsets of coupled perturbations with no dynamical 
mixing between the two subsets^. These are the following 

{^,e,/,a4 and {A„AX,b,B}. (27) 

Thus, we make no mistake ignoring the latter subset in this paper to study exclusively 
the first one. Furthermore, it was shown in [19] that the perturbations in the above gauge 
(which is called in [19], longitudinal gauge) coincide with the gauge invariant perturbations. 
Thus, there will be no gauge ambiguities in our results. 

•^The proof in [19] does not include the scalar modes, b and B, of the Kalb-Ramond field Bmn- However, 
it can be easily shown that for angle- independent perturbations, they have no mixing with the subset of 
V', / and ae. 
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3.1 Linearized equations of motion 



The linearized Ricci and energy-momentum tensors resulting from the above metric and 
field perturbations are given in Appendix A. Using them, we can write down the various 
components of the linearized Einstein equations. The {in/) component reads 



Dip 1 W W', la', q 

— + -lb" + 3 ih' + — f' + lb' + ^ 

2W^ 2^ 2a^ 



2/- ' 



8W^ 



/ = 



The (lip) component reads 
V 

The (pp) component reads 



W a' W 

^' + e' + 2— (^ + 0-2-(^-0-2-^/ 



q 



= . 



(28) 



(29) 



^ + ^b" + -f" + 4 — V + 2 — ^ - 2— V/ H £ - 2 f 

2W^^ 2^ a 2a W' 

, , 1 



8W^ 



The (^^) component reads 



— ^ - -0" + -£" - 4 — ib' + 2 — e + $ 

^ 2^ 2 a^ 



3g / 



2/ + 



3f_ 
8W^ 



, , 3g^ 1 



In addition, linearizing the scalar equation (5) we obtain 



3g 



2 I 3g^ 



3g2 



/ - . (30) 



/ = 0. (31) 



V[/2 ^-^ ^ a^ W^^ 



2 H 



2Vr8 



,2 H 



f- 



Q 



, (32) 



while from the gauge field equation (6) we get 

□a. 



— ;r + a/) + — an af) H 



f 

-3*' + I 



. 



From (29) we can solve for ae as a function of the other three perturbations 
XaW^ 



W a' W 

^'+r+2-(^+o-2-(^-o-2-^/ 



(33) 



(34) 



Since (33) does not have any singular (5-function) term, we should have at the bound- 
aries that a'g{ps) = 0. Then, from (34) we find that there exist two possibilities: 

(i) If ^ = then the only condition is that the derivatives of the wavefunctions should 
vanish at the poles, e.g., ip^Ps) = 0. 

(ii) If ^ 7^ ^, then we find the stronger condition that the wavefunctions themselves 
should vanish at the boundary with a limit for e.g., ip to be 

hm 44 = ■ (35) 

p^ps a[p) 

From the above, we see that since a{p) oc (p — p^) + . . ., the wavefunction in case (ii) should 
have an expansion around the singularities as'0 = C(p — Ps)^-|-..., which also implies that 

i^'iPs) = 0. 

Two of the remaining five equations are trivial. The first trivial combination is (pp) + 
{99) + 2(pz/). The second trivial equation is the gauge equation. To show the latter, we 
substitute ag from (34) to (33) and in the resulting expression we substitute □/ from (32), 
and D^', D^' from the derivatives of (28) and (30) respectively. When we do that, we find 
that some local terms (5-function terms) remain, namely 

[WW\f + i-2,i^)-,p']5{p-ps) (36) 

which give zero contribution since W'{ps) = and ip'{ps) = 0. 

Thus, finally we have three remaining equations for the three variables ip, C,, f ■ So, 
after we have substituted ag from (34) and simplified the expressions by means of the 
background equations of motion, we can rewrite Eqns. (28), (30), (32) as 

2^ 2^ 2a^ 2a W' 

M/2 2^ 2^ 2a^ 2a^ W-^ 

WW a a a 

[The (5-function terms in the above equations can be dropped for both cases of boundary 
conditions (i) and (ii).] 

In order to find the spectrum of the above system, we have to determine the relation 
between the perturbations which give rise to a single fluctuation equation. To do this, we 
can rewrite the system as 
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ox = rri^X , (40) 

with X — (V',^,/) and HX — m^X. Although the operator O does not look Hermitian 
in this X basis, we can consider the eigenvalue problem by looking for the linear relations 
between the perturbations which are consistent with the linearized equations. Substituting 
the following ansatze 

i{x,y) = A{y)tlj{x,y) and f{x,y) = B{y)i){x,y) , (41) 

to the above equations, we have to find the coefficients A and B that consistently collapse 
the system to a single differential equation. In this procedure, we find two distinct cases, 
in which both A, B are found to be constant. 

Constant wavefunctions 

The first case is the one where the wavefunctions are constant, i.e., ip = iIj{x). Then, 
it is easy to see from the above system, that there are two possible solutions for (^4, B) 

{A,B) = (1,2) and {A, B) = {1, -2) . (42) 

Thus, we have C = ip (therefore the singularity condition (i) is satisfied trivially) and in 
addition two possibilities for / = ±2-0. The first mode corresponds to the massless mode 
with Ht/j — and the second one to a massive mode with Ht/j — Idg'^t/j. For following use, 
we will call the first mode ■0o and the second one ipi. The above shows that the two modes 
which one finds in the unwarped case [17], maintain their form even when we introduce a 
warping. 

Let us note here that the massless mode, which corresponds to the breathing mode of 
the internal space, has the same relative wavefunction as the graviton zero mode. By this 
we mean that in the four dimensional part of the metric we have 

dsl = {W\p)rj^, + W\p)[h^,{x) - iP{x)rj^,]}dx^dx'' . (43) 

This is in contrast with the five dimensional case of e.g., the Randall-Sundrum model, where 
the relative wavefunctions of the radion and the zero mode graviton were different [20] . 

Non-constant wavefunctions 

The second case is the one where the wavefunctions ip{x,y) have non-trivial profiles. 
Then, one can see that the only way that all three equations (37), (38), (39) collapse to 
the same second order equation for is when 

(AS) = (-1,-2), (44) 

in other words, when ^ = —ip and / = —2ip. In this case we obtain the differential equation 
for the ffuctuation 
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As discussed before for the case (ii) of boundary conditions, the wavcfunctions at the 
poles of the internal manifold should satisfy ip'{ps) = '^{Ps) = and have an expansion 

^^;^c{p-p,r + .... 

Summarizing, the spectrum of the scalar excitations of the model consists of a zero 
mode, a first excited state with constant wavefunction and a tower of additional excited 
states with non-constant wavef unctions. In the unwarpcd case, the gauge field perturbation 
is zero for the zero mode and the massive constant wavefunction state, but non-trivial for 
all the other modes. On the other hand, in the warped case, the gauge field perturbation 
is nontrivial for all the massive states. 



3.2 Solutions for the non-constant modes 

To solve Eqn. (45) for the non-constant modes we will first separate variables as iplx, p) = 
"ip{x)x{p) with □■^ = m^%l). Then Eqn. (45) becomes 



Then, let us do the following transformation 



dp^Wdz , X=[-^,\ X. (47) 



1/2 



where we define 



a„.^sin(|). (48) 

Let us note here that the z and the p coordinates coincide for no warping {W = 1), and 
then additionally a = and X — X- With the new variables and wavefunctions substituted 
in Eqn. (45), we find that 

" do X /Id 3d^ 15 VF^ 1 dn 3dn\ . 9. „ , 

X--X+ o--7--XT^- o +7^ X + ^'X = , 49 
Oo \ 2 a Aa^ 4 2 oq 4 Oq / 

where ' = d/dz. Prom Eqns. (9), (11) and (17), we know the coordinates are related as 

r^rotan^^V (50) 



Then, using 

1/4 

/i ' " \h) 



a = r^ , W=(fy) , (51) 
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Ag/q ^ 1 Ag/q = 2 Ag/q = 1/2 




Figure 1: The wavefunctions x{p) for first three non-constant modes and for three different 
ratios of Ag/q. The first state wavefunctions are plotted with thick hues, the second with 
thin hues and the third with dashed hues. 



the terms in the parenthesis in Eqn. (49) cancel and we are left with the Schrodinger 
equation 

X X + m x = . (52) 

As seen from the above equation, in the ^-coordinates, the eigenvalues of the system depend 
only on Tq {i.e., on g) and are independent of the warping. The eigenfunctions x are as well 
independent of the warping, however, the wavefunctions % in (47) depend on the warping. 

The solution of this equation, with the boundary conditions previously mentioned, is 
given in terms of the Legendre's associated polynomials as: 



Xi 



^ 'l-y^f'^PlM with y^cos{^y^. (53) 



2l{l - 1) 

On the other hand, the eigenvalues of the equation are given as 

,2 4 



mf = —l{l-l) with Z = 2,3,4,- •• (54) 

Here, let us stress that the first level, for / = 2, has mass m2 = Ag, which is the same as 
the one of the constant massive mode ipi. Thus, these two degenerate states can mix, as 
we will see later. The variable y is related to the original r coordinate as 

1-4 

The Legendre's associated polynomials are given as 

- (1 - vr'-'^/M - ^^^^i^ - D" . (56) 
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satisfying the orthogonality condition 

dy P^(y)P„^(y) = ^-^^^ ■ . (57) 
Therefore, we obtain the following orthogonality condition between scalar modes, 

J ao J a 

The wavefunctions of the modes with Ag/q> 1 are localised closer to the brane sitting 
at r = oo, while the modes with 4:g/q < 1 are locahsed closer to the brane sitting at r = 0. 
Examples of these modes are given in Fig.l. 



4 Effective action for scalar modes 

Wc would like now to calculate the quadratic four dimensional effective action of the 
perturbations that we considered in the previous section. In order to do that, we have to 
expand the bulk and brane actions (1), (2) up to quadratic orders of scalar perturbations, 
substitute the modes that we have already found and integrate the extra two dimensions. 
This task is done in detail in Appendix B, leading to the following quadratic effective action 
for the scalar modes 



V'onV'o + ^V'i(n-i6r)V'i 



+B{i/j^Dil>2 + ^sD^/^i - 32^Vi^2) + CJ2M0- m^i)^i 

l>2 



l>2 



where Mj, = Xttv^ is the effective four dimensional Planck mass. The constants appearing 
in the action depend on the ratio 4:g/q and read 



A = 



B 



C = 



6 ^ 12 V g 



2v^L 
q ) 



1 + 



+ 



(60) 
(61) 
(62) 



In the above action, we also added the coupling of the scalar modes to the brane matter 
by the brane energy-momentum tensor T^l^"^^. 
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Since ipi and ■02 are degenerate in mass, i.e., — lOg^^, we can choose a basis such 
that there is no mixing term between these two states in the quadratic action. This is 
given by 



V2 



with 



d 



C-A 
2B 



Thus, we obtain the effective action in a diagonal form as 



(63) 



(64) 



eff 



2Mj. 



[^oD^o + i^+0+(n - 165^)0+ + K_ 



l>3 



(0+ + 0_)]<'^r; 



brane 



with 



(65) 



(66) 



VTTd^ 

These kinetic term coefficients are always positive definite. Therefore the solutions that 
we described are stable for any value of 4:g/q. 

Note also that there is no coupling of higher KK modes to the brane matter. So, from 
the brane perspective, we can think of the four dimensional effective theory, only consisting 
of the fields {ipQ, 0_). The two massive modes 0+ and although they are degenerate 
in mass, they have different couplings to brane fields, since when canonically normalised 
their couplings will get a contribution from the different and K^. 

In particular, for the unwarped solution, we have A = C = 1 and S = ^. Then, the 
effective action for the unwarped case is 



£eff = 2M|[z^on^o + (l + ^)0+(n-16/)0+ + (l--^)0-(n-16<72)0 



l>3 



/lI/^^brane 
'I -'■ ixu 1 

and the kinetic terms for the scalar modes are positive definite. 



(67) 
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5 Conclusions 



We have analyzed the hnearized scalar perturbations (which couple to the dilaton) of 
the general axisymmetric vacuum of the Salam-Sezgin model. We have found that the 
mass eigenstates consist of a zero mode with constant wavefunction, a degenerate pair 
of first excited states, one with constant and the other with non-constant wavef unctions, 
and a tower of heavier states with non-constant wavef unctions. The degenerate pair has 
quadratic mixing in the effective action, even in the case with no warping. The orthogonal 
combinations couple to the branes at the two poles with different strengths. It is important 
that their kinetic terms are positive definite for any value of 4,g/q and thus the model is 
stable. These results are expected to hold also for the anomaly free-model [3] as long as 
the additional scalars in the added hypcrmultiplcts arc much heavier. 

The above-mentioned extra degenerate mode that we have found has to be included 
in the effective four dimensional A'" = 1 supergravity description of the unwarped solution 
without branes [17]. The new chiral multiplet which has to be included will have a complex 
scalar component whose real part will be the extra degenerate mode. The imaginary part 
will be provided by the scalar perturbation sector that we did not study in this paper, most 
probably from the Kalb-Ramond field. A re-analysis of the effective theory is necessary to 
capture the influence of the new mode in the low energy physics. 

Furthermore, it would be interesting to study again the cosmological properties of this 
model, as in [21], including the non-constant mode. The analysis will be more complicated 
since the new mode dependence on p will complicate the derivation of a four dimensional 
potential. We plan to address this issue in a future work. 
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Note added 

In the original (and published [22]) version of the paper, there has been a typo in (63) 
and (66), which indicated an instability of the warped solutions for a certain region of the 
parameter space. This typo was noted in Appendix B of [23], and was corrected (along 
with the conclusion related to it) in the present manuscript. The relevant erratum has been 
as well published in [24]. Furthermore, it seems that our requirement that the system (40) 
collapses to a single differential equation is too strong and, thus, there exist modes [23] 
consistent with our boundary conditions which we have missed. The correct treatment of 
the system (40) reveals that the spectrum consists of a massless ground state and a tower 
of massive modes with triple degeneracy at each level. 
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Appendix A: Linearized equations for scalar modes 

The linearized Ricci components computed from the perturbed field ansatze (24), (25), 
(26) are 



HP 



Rpp 
Ree 



-W^ h3 — 

W aW 

1 

2 



W" W'^ a'W 

+ 3-77Tr + 



W 



aW 



(A.l) 
(A.2) 



-4 



W «" 1 ,// 1.// ,W W', a', 3a', 



a" 1 , 1 1 ,, ^Py' , , 3 a' , 

4^ TTT^nV^ + 7^T77?ne - V'" + of - 4— V^' + 2— + --r 



aiy 

II 



a 



21^2 



W 



W 



2 a 



(A.4) 



where □ = 'q^'^d^di,. The linearized energy-momentum tensor on the other hand is 



T, 



p,v 



T, 



HP 



T, 



pp 



Tee 



2W^ ' 27rAa' 



3q^ 
2W^ 



4W^ 



TlXa 

W q 

1 



2 





1 




2W^ 







1^2 
1 

1 



2 , 9 



4g^ + 



f + 2^f' + ^a' 



-4W"-8g^+ ^ 



21^8 



^ [8iy'2 + 16/] + ^ 



(A.5) 
(A.6) 



8W'^ + 8g^ - 



(A.7) 
2 n 



1^2 



2 , 9 



4/ + 



/-2t77/' + 



2W^ 
(A.8) 



The Einstein equations are in our conventions 



R 



MN 



Tmn — -^TgMN 



(A.9) 
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Appendix B: Quadratic action for scalar modes 



By expanding the bulk action (1) and the brane action (2) up to quadratic orders of 
scalar perturbations, we get the four dimensional effective Lagrangian as 



with 



■'int 



£eff ^2nJ dp\aW\Ck + + Ant) , 



fB.ll 



Ck = w- 



-v^r - e - 3^ - (V' + i^r + e + i2^V'' + + s^e' 



(B.2) 



1 „,o .VT' 



(-3V' + ^/)a', 



4^T^ + 4- 



2 27rAa 



(B.4) 
(B.5) 



We can check that the variation of this action reproduces indeed the correct linearized 
equations. 

Now, we can write the scalar perturbations as a sum of the modes we have found in 
section 3 



f 

ae 



n 

ipo {x) + -01 {x) - ^ (a;, p) , 

n 

2V>oW-2?),(i)-25^'i/.„(i,p) , 



= 4" + E' 



,(«) 



(B.6) 
(B.7) 
(B.8) 
(B.9) 
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with 



AXaW^ /W^ _ a' 



(1) _ o^avv vv („) 



a J 



(B.IO) 



Then, inserting the above perturbations into Eqn. (B.l), we can rewrite the 4D effective 
Lagrangian as 



nonconst ~r •'-mix ~r •'-int / ) 



with 



(B.ll) 



U = 2W- 



+2(^1 + —' 



W 

w^^°)(V'inV'n + V'nnV'i) 



8 VT' /FT' 



a 

i'n 2 



■-const 



.2 . 2 



'nonconst 



g2 V 1^ a / 



2?^ 



(B.12) 



(B.13) 



-2^/^i 



2 VW^ a 



VF' a' 

Ant = V"'(V'0 + V'l+V'nX"T, 



W^2 V 



hr Kp- Ps) 

27rAa 



[Here we note that the summation over indices for massive modes is understood.] 
Now let us simphfy the terms in the action. First, by using the fact that 



we obtain that the mass terms for the constant modes reduce to 



(B.15) 
(B.16) 



(B.17) 



■-const 



3V 



1 + 
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Moreover, using the fact that 



W a 



W a 



3V 



W a' 



(B.18) 



W 



10 



\W a) 2nXa 



S{p-Ps) , (B.19) 



integrating by parts the first and second terms in Eqn. (B.14) and using the equation of 
motion (46), we obtain 



'nonconst 



Am 



a / 

/W a'\- 
\W ~ "a/ 



10 



(B.20) 



The first term in Eqn. (B.15) vanishes as a surface term. So, after integrating by parts for 
the remaining double derivative term in Eqn. (B.15), we get 



By integrating the first term by parts, we finally obtain the mixing mass term as 



C. ■ = 



Q^W\W a) J 



(B.21) 



(B.22) 



Finally, making a separation of variables as ^„(x,p) = V'n(a;)Xn(p), we can rewrite the 
effective Lagrangian for scalar modes as 



'eff 



2M| 



V'onV'o + ^V'i(n-i6/)V'i 



ibr 



where Mp = Xiirl and 



A = 



Br, 



nl 



27r 

m| 

Ml 
Ml 



/''''--=(-?^(^-0-'")-' 

dp\aW^(l + 



2 _ a'\2 



(B.23) 



(B.24) 
(B.25) 
(B.26) 
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It can be seen through some algebra, that the factor multiplying XnXi (B.26) is the 
proper weight appearing in (58). Similarly, the factor multiplying Xn in (B.25) is the 
wavefunction of the lowest non-constant mode X2 times again the proper weight. Thus, 
from the obtained solutions for xi iii Eqns.(47) and (53), the above expressions can be 
simplified as following 
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